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Measurements of the nonlinear flux-flow resistivity p and the critical vortex velocity v* at high 
voltage bias close to the instability regime predicted by Larkin and Ovchinniko\"- are reported 
along the node and antinode directions of the d-wave order parameter in the a-b plane of epitaxial 
Y Ba2CuzOT-s films. In this pinning-free regime, p and v* are found to be anisotropic with values 
in the node direction larger on average by 10% than in the antinode direction. The anisotropy of p 
is almost independent of temperature and field. We attribute the observed results to the anisotropic 
quasiparticle distribution on the Fermi surface of Y Ba2Cm07-s. 
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Flux flow due to the Lorentz force in a current carry- 
ing type-II superconductor in the mixed state has been 
studied intensively for many years^ii. In view of the 
d2.2_j,2-wave symmetry of the pair wavefunction in the 
hole-doped cuprates^, measurements of the flux-flow re- 
sistance (FFR) as a function of the crystallographic di- 
rection within the Cu02 planes become highly interest- 
ing. In c-axis oriented epitaxial Y Ba2Cu30j^s (YBCO) 
fllms under an applied magnetic fleld parallel to the c- 
axis, one could envision an anisotropic behavior resulting 
in some difference between the flux motion in the node 
and antinode directions. However, up to now no such 
anisotropy has been reported, nor has such an effect been 
analyzed theoretically. Recently, a search study for flux 
flow resistivity anisotropy in the low bias pinning regime 
was carried out by some of us in epitaxial thin fllms of 
YBCO^. To within the experimental error in that study, 
no anisotropy was found between the transport proper- 
ties of microbridges patterned along the node and antin- 
ode directions. Since at low bias the transport results 
are determined by the pinning properties of the fllms, 
it was concluded that the pinning properties in YBCO 
are isotropic. At high bias however, close to the Larkin 
Ovchinikkov (LO) instability regime^, pinning is negligi- 
ble and flux flow resistivity measurements could reveal 
an intrinsic anisotropy between the node and antinode 
directions. In this study we present the first experimen- 
tal results indicating such anisotropy in the a-b plane of 
YBCO films, where the nonlinear FFR along the node 
direction is enhanced on average by about 10% as com- 
pared to that along the antinode direction. 

Our measurements were performed on the same two 
wafers used previously in the low bias measurements^. 
We used two high quality epitaxial, c-axis oriented 
YBCO films of 0.12 /im thickness prepared under iden- 
tical conditions by laser ablation deposition on (100) 
SrTiOs (STO) wafers of 10 x 10 mm^ area. The films 
were patterned with a photolithographic mask into 10 
equally spaced 0.12 x 12 x 100 fim^ bridges along the 
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FIG. 1: (Color online) A schematic drawing of the bridges on 
the wafers, and a node and antinode bridges voltage V versus 
current I at B = 6 T and T = 73 K. Inset (a) shows the 
corresponding I versus V just below the maximum of I versus 
V with fits using Eq. (1), inset (b) shows the dynamic nature 
of R\ and i?2, see text. 



line dividing the wafer into two halves by deep UV pho- 
tolithography and Ar ion milling. The dimensions of the 
resulting bridges were measured with optical and atomic 
force microscopes and found to conform with the mask 
design to within ±0.4%. Successive microbridges were 
oriented at alternating angles of 0° and 45° to the edge 
of the wafers as shown schematically in the upper panel 
of Fig. 1, so that the transport current would flow either 
along the node or the antinode directions of the order 
parameter. The alternating direction of adjacent bridges 
is important in order to minimize systematic differences 
due to possible large-scale inhomogeneities in the fllms. 
On one fllm, where the orientation of the STO substrate 
was with the edge of the wafer parallel to the (010) crys- 
talline direction, five odd number bridges were along the 
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antinode direction, and five even number bridges along 
the node direction. In the other film where the side of 
the wafer was parallel to the (110) orientation, the role 
of the antinode and node bridges was reversed due to the 
epitaxial growth of the film. Studying these two types 
of wafers was done in order to check if our ion milling 
process, done at an incident angle of 45° to the wafers, 
is affecting the properties of the bridges. Any observed 
difference in the transport properties of the two wafers 
would imply that the effect is not intrinsic, and results 
from the patterning process. Low resistance gold con- 
tacts were used, and the wafers were annealed for optimal 
oxygen doping of the films. The resistance measurements 
were done by the standard four-point dc technique in He- 
lium gas environment with and without a magnetic field 
of up to 8 T normal to the wafers. 

The low bias normal-state resistivity p under zero- 
magnetic-field was measured versus temperature for each 
bridge on the two wafers. Averaged values of the inner- 
most bridges are given in Figs. 1 and 2 of Ref. 5. The 
onset of the superconductive transition at B = was 
identical for the two kinds of bridges on the two wafers, 
with r°"*'=* = 94 K. This value of r°""'=* has also been 
reported by Gross et al& for the same kind of films. The 
resistivity in the normal state at 95 K was p{95 K) = 
168 fiflcm, being identical for both wafers and both type 
of bridges. The observation that p(node) = p(antinode) 
results from the four-fold symmetry and the heavy twin- 
ning of our films. 

At lower temperatures, the FFR in the pinning regime 
of YBCO was recently found to be isotropio^. To de- 
termine the pinning-independent FFR we measured the 
current-voltage characteristic near the LO instability 
caused by the nonequilibrium distribution of the quasi- 
particles at moderate electric fields^. This instability 
has been studied in detail both in conventional super- 
conductors and in the cuprates?. Following procedures 
of earlier experiments^, for each temperature and mag- 
netic field, the voltage applied to the sample was ramped 
up until a resistive jump due to the LO instability was 
reached. A typical such jump is seen in the TV curves of 
Fig. 1. Above this jump a thermal runaway is observed 
(not shown here), where the bridges reach the normal 
state with a typical resistivity of 200 /zJlcm which corre- 
sponds to about p{105K). We normally did not expose 
our bridges to this regime to avoid possible thermal dam- 
age. Eq. (53) of the LO theorjii yields 

a) 

where c is a number of order unity, V* is a critical 
voltage, and i?i and i?2 are two FFRs which were taken 
equal in the original theory. We stress that Eq. (1) is 
valid only in the flux flow regime at high bias near the 
maximum of I versus V, and not in the pinning regime 
at low bias. In the low bias regime where / is linear in 
V, the FFR was found to be isotropic^, as is expected 



for a heavily twinned YBCO crystal with a four-fold 
symmetry if the pinning is also isotropic. By fitting our 
I-V data in the instability regime using Eq. (1) we find 
that the FFR Ri/{1 + y/l - T/Tc) and R2 have similar 
values to within 10%. Since the Lorentz force is oriented 
perpendicular to the current and the applied magnetic 
field, and YBCO has a four-fold symmetry, flux flow is 
along the node direction in the node bridges, and in the 
antinode direction in the antinode bridges. The resistive 
voltage across the bridge is F = x B|L, where v^p is 
the flux-flow velocity, B is the magnetic fleld and L is 
the bridge length. The critical voltage V* in Eq. (1) is 
thus associated with a critical vortex velocity v* by this 
relation. 

Joule heating of our samples could affect our results. 
The possible temperature rise is AT = aP, where P is 
the dissipated electric power per unit area of the sample, 
and a is the thermal boundary resistance®. YBCO fllms 
on SrTiOs substrates have a = lO^-^cm^K/W at 300 
K®. Assuming an even larger value of 10^'^cm^K/W, 
and taking the maximum values of P close to the LO 
instability, we estimate that AT remains smaller than a 
few tenth of a Kelvin. We therefore conclude that Joule 
heating of our samples can be neglected, except perhaps 
close to Tc. Experimentally, different rates of the current 
sweeps (14 - 300 mA/s) showed reproducibility of the 
I-V data, thus excluding the possibility of magnetic 
hysteresis^ and thermal heating effects. 

Two typical I-V curves for a node and antinode 
bridges are shown in Fig. 1 at T = 73 K and B = 6 T. 
The difference between the two curves is very clear and 
typical of most of our bridges. With increasing current 
above the critical current regime, the curve becomes 
steeper and steeper up to the instability. In the case 
of voltage-controlled operation as we have, a negative 
differential resistance sets in just before the jump as seen 
in Fig. 1. In the pinning-independent instability regime, 
the flux-flow process follows closely Eq. (1). Indeed, in 
inset (a) to Fig. 1 where ten data points of / versus V 
just below the maximum of / are plotted for each curve, 
one can see that the fits using Eq. (1) are in excellent 
agreement with the data. In the fits we assumed that 
c = 1, and obtained R2 and V*. In the following we 
use the procedure as in inset (a) to Fig. 1, to present 

our results. 

Plots of the flux flow resistivities pi/{l + \/l - T/Tc) 
and p2 (calculated from Ri and R2) versus the magnetic 
field B are shown in Figs. 2 and 3 for the two wafers, 
respectively, at various temperatures. The spread of 
the mid-point Tc (where the low bias dR/dT is highest) 
of the different bridges on the two wafers are given in 
the insets to these two Figures. The missing data for 
bridges #4 and #9 in Fig. 2 is due to bad contacts. The 
relatively large spread of the Tc values on this wafer of 
about 1.2 K, and the asymmetry relative to the center of 
the wafers is due to an off-center alignment of the laser 
ablated plume on the wafer. Nevertheless, bridges #5-8 
on LAF551 (Fig. 2) and bridges #3-6 on LAF552 (inset 
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FIG. 2: (Color online) Averaged flux flow resistivity pi/(l + 
i/l - T/Tc) versus magnetic field B on LAF551 at 60, 70 and 
78 K, for the node (stars) and antinode (squares) directions. 
The inset shows the distribution of the Tc values on the wafer. 



FIG. 4: (Color online) Histogram of the anisotropy results on 
both wafers Ap2=p2(node)-p2 (antinode) normalized by the 
average value of p2 for all fields and temperatures. 
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FIG. 3: (Color online) Averaged flux flow resistivity p2 versus 
magnetic field B on LAF552 at 67, 73 and 77 K, for the node 
(stars) and antinode (squares) directions. The inset shows the 
distribution of Tc for the four measured bridges on the wafer. 

of Fig. 3, in a reversed order of the bridges) have a much 
smaller spread of about 0.4K. Thus the FFR in these 
bridges was measured and compared, and no systematic 
correlation with the corresponding Tc values was found. 
In Figs. 2 and 3 average values are shown of the two 
measured node bridges and the two measured antinode 
bridges on each wafer. The errors in these figures are due 
mainly to this averaging procedure. These figures show 
that /9i/(l + ^/X^^^Tflc) and P2 increase almost linearly 
with B at low temperatures, but more than linearly at 
higher temperatures. In all cases, a clear anisotropy is 
observed where the pi/(l + ^1 — T/Tc) and p2 values 
along the node direction are larger than along the 
antinode direction. The fact that the linear term in 
Eq. (1) i?i oc p\ was found to be anisotropic seems 



to contradict the four-fold symmetry of the twinned 
crystal. We note, however, that R\ (and also i?2 and 
V*) are voltage dependent (see inset (b) to Fig. 1). The 
results of this inset were obtained by fitting the antinode 
I-V data using Eq. (1) near the instability regime at 
various voltage segments of ± 5 mV. The segments 
contained between 6 data points (above the maximum 
of I and just below the voltage jump) to 17 points at the 
lowest voltage bias shown. All fits had of better than 
0.999 (solid symbols). Also shown, are results of two 
inferior fits with 48 and 65 data points (empty symbols), 
and as one can see, the resulting parameters still follow 
the general behavior. We thus conclude, that i?i is not 
a linear term of the FFR, and therefore does not have to 
obey the symmetry of the heavily twinned crystal. 

Fig. 4 shows a histogram of the normalized 
anisotropy A/92=P2(node)-/92 (antinode) for all the p2 
data on the two wafers. In this figure Ap2 is normalized 
by the average value [p2 (node) -I-/92 (antinode)] /2 at each 
field and temperature. We first plotted the normalized 
anisotropy Ap2 versus field for the six different temper- 
atures, and found that it is basically field independent 
to within the noise of the measurements. No systematic 
temperature dependence of the anisotropy was found at 
60, 67 and 77K as can be seen in the histogram of Fig. 4. 
This apparently is due to noise in these measurements. 
At 70, 73 and 78K however, one can see that the 
anisotropy values are concentrated in a narrow range of 
about 11 ± 2%. We thus conclude that the anisotropy 
is temperature independent in the temperature range 
of our measurements. We note that the anisotropy 
is expected to vanish at high temperatures due to 
increasing thermalization of the quasiparticles. Above 
80 K however, we could not observe the LO instability 
jump in the I-V curves, to observe this effect. 

Next, we turn to the critical voltage V* which 
corresponds to the critical vortex velocity v* according 
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FIG. 5: (Color online) Critical vortex velocity v^ calculated 
from V* of the fits to Eq. (1) versus B for the node (stars) 
and antinode (squares) directions for the bridges of Fig. 2. 
The lines are guide to the eyes. The inset shows a schematic 
drawing of the twin boundaries in the two kind of bridges. 

to V* = |v* X B|L. In Fig. 5 we plot average values of 
the critical vortex velocity v* of the node (#6,8) and 
antinode (#5,7) bridges on LAF551 versus magnetic 
field. In the node direction is systematically larger 
(by about 10-20%) than in the antinode direction. In 
LAF552 this anisotropy is similar (not shown). The 
V* values at each temperature show only a weak field 
dependence, and are basically quite constant as was 
found previousljii. Our conclusion therefore is that the 
anisotropy in is consistent with the anisotropy in pi 
and p2, as all have comparable percentage values. 

In the absence of a theoretical prediction for the 
anisotropics of the flux flow resistivities Api and Ap2, 
and the vortex velocity Av*, we shall present simple 
qualitative arguments for explaining them. In principle, 
in the absence of pinning, the vortex velocity v* is 
affected by the ability of the vortices to transfer momen- 
tum to the quasiparticles. The Magnus or Lorentz force 
that acts on a vortex constantly transfers momentum 
to it. After a short acceleration, the vortex reaches a 
terminal velocity v* in the viscous media while the extra 
momentum has to be constantly dissipated to the quasi- 
particle excitations. Thus the easier it is for a vortex to 



transfer momentum when it moves in a certain direction, 
the faster it will move in this direction. It was predicted 
theoretically that in a d-wave superconductor the low 
energy excitations are located in the vicinity of the 
nodes'^'^ . This was verified in photocmission experiments 
where four small regions of quasiparticle excitations 
were found on the Fermi surface in the directions of 
the nodesii. Our experimental observations indicate 
larger flux flow velocity and FFR in the node direction. 
Therefore, it seems that the momentum transfer by 
vortices to the quasiparticles is more effective in the 
node direction as compared to that along the antinode 
direction where less quasiparticles are available for the 
momentum transfer process. Clearly, a comprehensive 
theoretical analysis of this issue is needed in order 
to make a detailed quantitative comparison with the 
present results. 

Finally, we rule out the possibility that twinning 
is responsible for the observed anisotropy. Twins in 
thin YBCO fllms form a dense mosaic of elongated 
crystallites with boundaries along the (110) and (110) 
node directions as shown schematically in the inset to 
Fig. 5. If the twin boundaries serve as easy channels for 
flux flow, then flow in the anti-node bridge is unhindered 
since there is always a Lorentz force component in the 
direction of the motion. In the node bridge however, the 
flow is blocked each time it reaches a boundary normal 
to its direction. This should lead to a larger flux flow 
velocity and FFR in the antinode direction as compared 
to the node direction. Since this is opposite to the 
observed results, twins could not be the source of our 
observations. 

In summary, we observed a clear anisotropy in the 
nonlinear flux flow resistivities pi and p2 of YBCO 
which we attribute to the intrinsic anisotropy of the low 
energy quasiparticle distribution on the Fermi surface 
of a d-wave superconductor. The observed anisotropy 
in the critical vortex velocity v* is consistent with the 
anisotropy of the measured flux-flow resistances. 
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